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Abstract 

\ We study a class of nonlinear pricing models which involves the feedback effect 

PLh ' from the dynamic hedging strategies on the price of asset introduced by Sircar 

. and Papanicolaou. We are first to study the case of a nonlinear demand function 

Lf^ \ involved in the model. Using a Lie group analysis we investigate the symmetry 

qh' properties of these nonlinear diffusion equations. We provide the optimal systems 

of subalgebras and the complete set of non-equivalent reductions of studied PDEs 
to ODEs. In most cases we obtain families of exact solutions or derive particular 
J> ' solutions to the equations. 
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1 Introduction 



X 

$H ' One of the important assumptions of the classical Black-Scholes theory is the 

^ assumptions that any trading strategy of any trader on the market do not affect 

asset prices. This assumption is failed in the presence of large traders whose 
orders involve a significant part of the available shares. Their trading strategy 
has a strong feedback effect on the price of the asset, and from there back onto 
the price of derivative products. The continuously increasing volumes of financial 
markets as well as a significant amount of large traders acting on these markets 
force us to develop and to study new option pricing models. 
There are a number of suggestions on how to incorporate in a mathematical 
model the feedback effects which correspond to different types of frictions on the 
market like illiquidity or transaction costs. Most financial market models are 
characterized by nonlinear partial differential equations (PDEs) of the parabolic 
type. They contain usually a small perturbation parameter p which vanishes if 
the feedback effect is removed. If p tends to zero then the corresponding nonlinear 
PDE tends to the Black-Scholes equation. 
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+ TTTiTT, ^=0' (1-1) 



Some of the option pricing models in illiquid markets possess complicated 
analytical and algebraic structures which are singular perturbed. Wc deal with 
singular perturbed PDEs if one of the nonlinear terms in the studied equation 
incorporates the highest derivative multiplied by the small parameter p. It is a 
demanding task to study such models. Solutions to a singular perturbed equation 
may blow up in the case /? = and may not have any pendants in the linear case. 
An example of a singular perturbed model is the continuous-time model developed 
by Frey [8] . He derived a PDE for perfect replication trading strategies and option 
pricing for the large traders. An option price u{S, t) in this case is a solution to 
the nonlinear PDE 

2 (1 - p\{S)SussY 

where t is time, S denotes the price and a the volatility of the underlying as- 
set. The continuous function \{S) included in the adjusted diffusion coefficient 
depends on the payoff of the derivative product. The Lie group analysis and 
properties of the invariant solutions to Eq. (1.1) for different types of the func- 
tion A(>S') were studied in [4]- [6]. The analytic form of the invariant solutions to 
this model allow us to follow up the behavior of these solutions. 

Under the similar assumptions Cctin, Jarrow and Protter [7] developed a 
model which includes liquidity risk for a large trader. Liquidity risk is the addi- 
tional risk due to the timing and size of trade. The value u{S, t) of a self financing 
trading strategy for the large trader in this setting is a solution of the following 
nonlinear PDE 

ut + ^a''S''uss{l-pSuss? = Q. (1.2) 

This equation seems to be simpler as the previous one but it is still a singular 
perturbed PDE. The Lie group analysis and the symmetry algebra admitted by 
this equation were studied in [2]. 

Sircar and Papanicolaou in [14] present a class of nonlinear pricing models 
that account for the feedback effect from the dynamic hedging strategies on the 
price of asset using the idea of a demand function of the reference traders relative 
to the supply. They obtain a nonlinear PDE of the following type 
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U-\l-pus)U'{U-\l-pus)) 



o rr-ui \TTl^TT-l^^ \\ q a S uss + r{Sus - u) = 

2 \ U ^{1 - pus)U'{U ^{1 - pus)) - pSussj 

(1.3) 

where t is time, S and a is the price and the volatility of the underlying as- 
set respectively, and the parameter r is the risk-free interest rate. The value 
u{S, t) is the price of the derivative security and depends on the form of the de- 
mand function U{-). The expression U^^{-) denotes the correspondingly inverse 
function, because of the strong monotonicity of the demand function the exis- 
tence of the inverse function U~^{-) is guaranteed. In the bulk of their paper 
[14] authors studied the particular model arising from taking [/(■) as linear, i.e. 
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U{z) = Pz, P > Q. The authors mainly focused on the numerical solution and 
discuss the difference to the classical Black-Scholes option pricing theory. 

In the present paper we study a more general case in which the demand 
function of the type U{z) = fiz°', a, /3 7^ is incorporated. Consistency of (1.3) 
with the Black-Scholes model characterizes the class of the admitted demand 
functions and leads to the condition U'{z) = Paz°'~^ > 0. In this case the model 
(1.3) takes the form 



du 1 
dt^2 



1-p 



du 
dS 
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The diffusion coefficient in Eq. (1.4) depends on both, us and uss multiplied by 
the small perturbation parameter p. It depends also on the parameter a which 
characterizes the type of the demand function and on the interest rate r. 

We give a short overview of analytical properties of Eq. (1.4) in the next 
Section 2. In Section 3 we provide the Lie group analysis of this equation. De- 
pending on whether the interest rate r = or r 7^ 0, we obtain different Lie 
algebras admitted by the respectively equation. Then wc provide optimal sys- 
tems of subalgcbras in the both cases. The optimal systems of subalgebras give 
us the possibility to describe the set of independent reductions of these nonlinear 
PDEs to different ordinary differential equations (ODEs) . In some cases we found 
the explicit solutions to these equations. We discuss the properties of invariant 
solutions in Section 4. 



2 Basic analytical properties of Eq. (1.4) 

In the model (1.4) introduced by Sircar and Papanicolaou in [14] the diffusion 
coefficient has a very complicated analytical and non trivial algebraic structure. 
In particular the diffusion term is represented by a fraction which contains deriva- 
tives us{S,t),uss{S,t). The authors analyzed some analytical properties of this 
equation in the case a = 1 in the vicinity of the Black-Scholes equation and con- 
sidered a valuation of a European derivative security with a convex payoff using 
their model. They give asymptotic results for when the volume of assets traded 
by the large traders is small compared to the total number of unites of the asset. 
In this Section we discuss some global analytical properties of Eq. (1.4) for a ^ 1. 

We pay the main attention to the second term in (1.4). Wc assume that 
the space variable 5 € U {0}, where ^2 = and the time variable t lie in 
T U {0}, where T = M"^. This term can vanish for some values of the variable S 
or on some set of smooth functions and then the equation may change the type 
from the parabolic one to another one. Other hand the fraction in the second 
term may became meaningless because of vanishing of the denominator on some 
set of smooth functions. We should exclude such functions from the domain of 
definition of our model. 

The classical linear diffusion equation of type Ut = uss is well defined on the 
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space D = C'^''^{n x r)nC({^ U {0}} x {TU {0}}) and u{S,t) map the space 
D to a space of continuous functions M = C({ri U {0}} x {TU {0}}). 

Let us check whether the expression for the diffusion coefficient in Eq. (1.4) 
vanishes or has singularities. 

Fist we study the case that the denominator of the fraction in (1.4) is equal to 
zero, i.e. we have to solve the equation 

1 -pus- -Suss = 0. (2.1) 
a 

It is easy to see that this equation has the following solution 

Using[S,t) = - + Ci{t) -S « +C2{t), a ^ 1, 

p a — 1 

Using{S,t) = ^ + ci(t)ln(5) + C2(t), a = 1, (2.2) 

where ci(t) and C2(t) arc arbitrary functions of t. We can rewrite the expressions 
in (2.2) as one expression which includes the case a = 1 as a limit case. Then we 
obtain 

Using{S,t) = -+Ci{t)-^ (S^ - l) +C2{t). (2.3) 

p a — 1 V / 

The numerator of the second term in (1.4) is equal to zero if one of the equations 
is satisfied 

S'^uss = 0, 

l-pus = 0. (2.4) 

The first equation is satisfied on all linear functions of S and in the point S = 0, 
the second equation has the following solution 

uoiS,t) = - + C2(,t). (2.5) 
P 

We notice that in the case ci{t) = the functions Using{S, t) and uo{S, t) coincide. 
It means in this case the numerator and the denominator of the fraction in the 
equation (1.4) are simultaneously equal to zero. 

In the second step we should define a limiting procedure to explain what we 
means if we say that (2.5) is a solution to (1.4). 

We chose in the space D a one-parametric family of functions u^ {S, t) of the 
following type 

u,{S,t) = di{t)S + d2{t) + ev{S,t), (2.6) 

where e G M is a parameter, the functions di{t),d2{t) are arbitrary functions of 
time and v{S,t) G D. If now the parameter e — > then the family of functions 
of the type (2.6) converges in the norm of the space D to a linear function of S, 
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i.e. to uo{S,t) = di{t)S + d2{t). We apply to this family the differential operator 
defined by (1.4) 



d[{t)S + d',{t) + ev,iS, t) + ^ eS-vssil - Pd.it) - epvs? ^ ^ ^^.7) 

2 (1 - pdi{t) - epvs - ePpSvss) 

here d-^{t),d2{t) denotes the first derivatives of the corresponding functions. From 
(2.7) follows that any linear function of S with constant coefficients will be solu- 
tion to equation (1.4) if the last term in (2.7) is a bounded function in the norm 
of the space M. If we replace the linear part in (2.7) by uq{S, t) we see that this 
function (2.5) is a solution to (1.4) if C2{t) = const.. The fraction in (2.7) is not 
well defined just in one case if v{S, t) coincide with the second term in (2.3). So 
far we use as the domain of definition for our model the space D the functions 
of type (2.3) with ci(t) ^ are excluded because they or their derivatives have 
singularities in the point 5 = and consequently they do not belong to the space 
D. 

In the classical case of a linear parabolic diffusion equation solutions of type 
(2.3) which do not belongs to the set of classical solutions are called viscosity 
solutions [1], [9] and these solutions arc well studied. 

We proved that the functions of type (2.3) should be excluded from the fur- 
ther investigation because the model (1.4) is not well defined on them. The linear 
function (2.5) with C2(t) = const, is a solution to (1.4) because any one paramet- 
ric family of functions Ue{S,t) in the norm of the space D convergent to uq is 
mapped by the differential operator defined by (1.4) to a zero-convergent family 
of functions in the norm of the space M. 



3 Symmetry properties of the model 

We provide in this Section the Lie group analysis of Eq. (1.4) first for the case 
r = then for r 7^ 0. In both cases it is possible to find the non-trivial Lie 
algebras admitted by the equation. Wc use the standard method to obtain the 
symmetry group suggested by Sophus Lie and developed further in [12], [11] and 
[10]. In the case r = we obtain a four dimensional Lie algebra L4 and by r 7^ 
Eq. (1.4) admits a three dimensional algebra L3 defined in the subsection 3.2. 

All three and four dimensional real Lie algebras and their subalgebras were 
classified by Pattcra and Winternitzs in [13]. The authors looked for classifi- 
cations of the subalgebras into equivalence classes under their group of inner 
automorphisms. They used also the idea of normalization which guarantees that 
the constructed optimal system of subalgebras is unique up to the isomorphisms. 

The symmetry group G4 related to the symmetry algebra L4 is generated by a 
usual exponential map. We use the similar procedure to obtain to each subalgebra 
hi from the optimal system of subalgebras the correspondingly subgroup Hi. 

The optimal system of subalgebras allows us to divide the invariant solu- 
tions into non-intersecting equivalence classes. In this way it is possible to find 
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the complete set of essential different invariant solutions to the equation under 
consideration. 

Using the invariants of these subgroups we reduce the studied PDE to different 
ODEs. Solutions to these ODEs give us the invariant solutions to the nonlinear 
PDE (1.4) in an analytical form. In the both cases whether by r = or by r 7^ 
we skip the study of invariant reductions to the two- and three- dimensional 
subgroups because of they give trivial results for the studied equation. 



3.1 Symmetry reductions in the case r = 

In the first step we solve the Lie determining equations for the equation 

^ (1 - pus - ^Suss) 

and obtain the Lie algebra admitted by this equations. We formulate the results 
in the following theorem. 

Theorem 3.1. Eq. (3.1) admits a four dimensional Lie algebra L4 with the 
following infinitesimal generators 

S d fS \ d d d ad d 

The commutator relations are 

[ei, 63] = [61,64] = [62, 63] = [62, 64] = [63, 64] = 0, [61,62] = 62, (3.3) 

Remark 3.1. In the very short letter [3] there is a misprint in the theorem 
formulation. We apologize by readers for the inconvenience. 

The Lie algebra L4 has a two-dimensional subalgebra L2 =< 61,62 > spanned 
by the generators 61,62- The algebra L4 is a decomposable Lie algebra and can 
be represented as a semi-direct sum L4 = i^2 © 63 ^ 64. The optimal system of 
subalgebras for L4 were provided in [13] and presented in Table 1. 

The optimal system of the one-dimensional subalgebras involves four subal- 
gebras h'^, i = 1, . . . ,4. We take step-by-stcp each of these subalgebras and 
the corresponding symmetry subgroup Hf and study which invariant reductions 
of the studied PDE are possible. 

Case Hf. This one-dimensional subgroup C G4 is generated by the 
subalgebra =< 62 > = < >■ It means that we deal with a subgroup of 
translations in the u - direction. Hence, to each solution to Eq. (3.1) we can 
add an arbitrary constant without destroying the property of the function to be 
solution. This subgroup does not provide any reduction. 
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1 


hi =< 62 >, /l2 =< cos (0) + 64 sin ((/)) >, 

hs =< ei + a;(e3 cos {(p) + 64 sin (0)) >, 
/i4 =< 62 + e(e3 cos {(p) + 64 sin (0)) > 


2 


/is =< 61 + x (63 COS ((?!)) + 64 sin (0)), 62 >, he =< 63,64 >, 
/ly =< 61 + a;(63 COS {(/)) + 64 sin (<^)), 63 sin {(/)) — 64 cos {(/)) >, 
/ig =< 62 + e(e3 cos (0) + 64 sin (</))), 63 sin {(/)) — 64 cos {(p) >, 
hg =< 62, 63 sin {(p) — 64 cos {(p) > 


3 


hio =< 61,63,64 >, /ill =< 62,63,64 >, 
hi2 =< 61 + a;(e3 cos {(p) + 64 sin ((?!>)), 63 sin (^) — 64 cos (<?!'), 62 > 



Table 1: [13] The optimal system of subalgebras hi of the algebra L4 where 
xeR, e = ±1, (pe [0,tt]. 



Case • The subalgcbra is spanned by the generator e3Cos{(p)+e4sin{(p). 
In terms of the variables 5', t, u it takes the form 

h^ =< ^ cos(0) + (^pS-^ + S-^^ sin(0) > . (3.4) 

The invariants z,w oi the corresponding subgroup H2 are equal to 

^; = S'exp(— fptan(<^)), w = u—-S (3.5) 

and we take them as the new dependent and independent variables, respectively. 
Then the PDE (3.1) is reduced to the ordinary differential equation of the fol- 
lowing form 

-p6zw, + la^w,,z^ ( ^ ] =0, S = tan{(P),^£[0,Tv],(P^^. (3.6) 

2 \ awz + zwzz J 2 

This second order ODE is reduced to the first order equation by the substitution 
Wz = v{z) which takes the form 

,.(_,i + !^,(l„„),(_^)')=0. (3.7) 

Eq. (3.7) has two trivial solutions z = and v{z) = Wz = which are not very 
interesting for applications. The non trivial solutions we obtain if we set the last 
factor in Eq. (3.7) equal to zero. We obtain the solution to (3.7) in the form 

^2 

wiz) = ciz^, p = l + a — ai: sj a(a — 2a), a = - — -, (3-8) 

Apd 

where ci is an arbitrary constant. In terms of the variables S,t,u the solution 
(3.8) is equivalent to the following solution to Eq.(3.1) 

^2 

u(S,t) = ciS^ ex.-p(—pptein((p) t), p = l + a — a± -v/ a{a — 2a), a = - — -, (3.9) 

4pd 
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where S = tan(^) , G [0, tt] , ^ 7^ § . 

Case H^. The subalgebra /13 is spanned by 

hi =< xcos{(p)^ + (xpsm{(j)) (^(^■^+xsm{(j))^ S-uj -^> . 

The invariants z, w of the corresponding subgroup are given by the expressions 

z = Se-'*, w = S^u- -S^+\ (3.10) 

P 

where b = {xpsm{(p) — ^)~^, c = (6xcos(^))~^, x 7^ 0, ^ € [0,7r],^ ^ ^ and 
xpsin(0) — i 7^ 0. 

We use z,w as the new invariant variables and reduce Eq.(3.1) to the ODE 



cr^a^ (z'^Wzz - '^bzwz + b{l + b)w) {zwz: - bwf 

-czwz H 2 ~ 0- (^-ll) 

2 (a {zwz — bw) + {z'^Wzz — 'ibzwz + 6(1 + 6) w)) 



Eq. (3.11) admits a solution of the type 

w{z) = cizi, (3.12) 

where g is a real root of the polynomial of the degree 5 

-cq{q{q - 1) + (a - 26) g + 6 (1 + 6 - a)f 
2 2 

+^{q{q - 1) - 269 + 6 (1 + 6)) {q - hf = 0. 



Eq. (3.11) has a complicate structure and is hardly possible to solve it in the 
general form. But by a special values of involved parameters we can simplify the 
equation and obtain some particular classes of solutions. 
We take the special case of Eq. (3.1) with a = 1. 

Under the special choice of the parameters (/> = 0, tt and 6 = —2 in Eq. (3.10) 
and by using the invariants z,w in the form 

z = Sexp ,xy^O, w = S~'^u - {pS)~^ , 

we reduce Eq.(3.1) to the ODE 

zwz{{z{z''w)z)zf + (T^x {z^w)zz{{z^w)zf = 0. (3.13) 

The substitution (3.12) in this case leads to the second order algebraic equation 
on the value of the parameter q 

q'^ + q{2 + a^x cos(^)) + a'^x cos = 0, q^ -2, ^ = 0, tt, (3.14) 
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which has two roots qi = —a'^xcos{(j)) and q2 = —2. For the future study we take 
just the first value qi. The value q2 = —2 leads to the know solution uo{S, t) (2.5). 
Since by q2 = —2 both, the numerator and the denominator in the fraction in 
(3.1) vanish, and the solution (3.12) coincide then with uo{S,t) by C2{t) = const. 
which we discussed in the previous section. 

We notice that the solution (3.12) differs from the function Using{S,t) (2.2) 
which involves the logarithmic term in the case a = 1. 

The solution to (3.13) or respectively to Eq.(3.1) in the form (3.12) in terms 
of S,t,u variables is equal to 

u{S, t) = ^ + Ci^^-'^'^e-'^* + C2, -2, (3.15) 

where Ci, C2 are arbitrary constants, a; G M and the first term is the only term 
which contains the dependency on the parameter p. We skip the factor cos(^) by 
X in this expression because cos(ij!)) = ±1 in the case = 0,7r and a; G M. 

It is remarkable that the reduced equation (3.13) does not contain any more 
the parameter p. Hence all invariant solutions of this class can be represented as 
a sum of two terms: the first one is equal to S/ p and the second one depends on 
z only but not on the parameter p. 

If we left the values of parameters like in the previous case, but take the 
invariants in another form z = InS* + tb/4: and w{z) = {u/S — l/p)S'^ than we 
obtain the different form of the rcdTiccd ODE 

w^^^w, + u;^^(u;2(4(l - 7) + k) + u;^u;(2(l - 7)^ + 2k(1 - 7)) 
+w^k(1 - jf) + wl{4{l - 7)2 + k{1 - 27)) + wlw{l - 7)(4(1 - 
+k{2 - 57)) + W;,w'^{l - 7)2((1 - 7)2 + k{1 - 47)) - K7(l - 7)W = q, (3.16) 

where k = 2a'^ /b. This second order equation can be reduced in the case Wz ^ 
to a first order ODE. We substitute p(w) = Wz{z{w)) and correspondingly Wzz = 
Pwp, i.e., w is the independent variable and p is the dependent variable in this 
case. Then we obtain the first order ODE 

pIp"^ + PwPipH'^i'^ - 7) + /t) +pw{2{l - 7)2 + 2k{1 - 7)) + w\{l - 7)2) 
+/(4(1 - 7)^ + '«(! - 27)) + P'^w{l - 7)(4(1 - 7)2 + k{2 - 57)) 

+pu;2(l - 7)2((1 - 7)2 + k{1 - 47)) - K7(l - jfw^ = 0.(3.17) 

This equation is quadratic in the first derivative Pu, and it is equivalent to the 
system of two first order equations. For some values of the constants 7 and k it 
can be explicitly solved. The simplest case we obtain if we chose 7 = 1 then the 
solution coincide with (3.15). In other cases the equation can be studied using 
qualitative methods. 

Case H2. We consider subalgebra h^^ spanned by 

hi =< €Cos{(p)-^ + epSsm{<p)^ + (1 + e5sin((^))^ > . 
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The invariants w are given by the expressions 



71" 

z = S'exp(-tptan(^)), ^ G (0,7r),^ 7^ -, 

w = — . , , , In 5* H — S — V, e = ±1. 

psm[(p) p 

Using these expressions as the independent and dependent variables we reduce 
the original equation to the ODE 



-ptSin{4>)wzZ + ^cr'^iwzzz'^ + a) 



a{wzZ — a) 



a{wzZ — a) + a + WzzZ^ 



0, (3.18) 



where a = (e/9sin(^)) ^ and cf) G (0,7r), {(/)) ^ |. Eq. (3.18) is possible to reduce 
to the first order ODE 



{v + a) {zvz — (1 — cx)v + a) 



-V [zVz — V + a) = 



2ptan{(f>) 



(3.19) 



after substitution v{z) = zwz — a. It is a quadratic algebraic equation on the 

value ZVz which roots depends on v only. If we denote the roots as f±{v) we 
represent the solutions to Eq. (3.19) in the parametric form 



dv 



Inz + ci, ci G M. 



(3.20) 



3.2 Symmetry reductions in the case r 7^ 

Theorem 3.2. Eq. (1.4) admits a three dimensional Lie algebra L3 with the 
following operators 



ei = 7^, e2 = (5-p)— , €3 = 8— +u-^ 



dS du' 



(3.21) 



The algebra L3 is abelian. 



Similar to the previous investigation for r = first we find an optimal systems 
of subalgebras for the algebra L3. Because we are interested just in the one- 
dimensional subalgebras we provide the optimal system for these subalgebras 
only 

hi =< 63 >, /i2 =< 62 + xes >, /i3 =< ei + xe2 + ye^ > . (3.22) 

First we provide for each of these three one-dimensional subalgebras hj. i = 1, 2, 3 
and the corresponding subgroup Hi, i = 1,2,3 a set of invariants. Then we use 
the invariants as the new independent and independent variables and reduce Eq. 
(1.4) to some ODEs. In the cases where it is possible we solve the ODEs. 
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Case Hi. The algebra hi is spanned by 

It describes scaling symmetry of Eq. (1.4) and means that if we multiply the 
variable S and u with one and the same non vanishing constant the equation will 
be unaltered. Respectively the invariants of this transformations are 

z = t, w = ^. (3.24) 

If we use these expressions as the new invariant variables we obtain a rather 
simple reduction of the original equation 

W;, = 0, (3.25) 

with the trivial solution w = ci = const. It describes all solutions to Eq. (1.4) of 
the type 

u{S, t) = Sw{z) = Sw{t) = Sci. (3.26) 

Case H2. The second subalgebra ^2 in the optimal system of subalgebras 
(3.22) is given by 

h2=<x S— + iS + u{x-p))—>, xeR. (3.27) 
ob ou 

The invariants of the subgroup H2 are 

z = t, w = + S^-^u, x^O. (3.28) 



If we use the invariants (3.28) as the new dependent and independent variables 
then the reduced Eq. (1.4) takes the form 

\ 2{ax — pY X J 

It has the solution 

w{z) = ci exp (p ( ^^^ + -1 V> ci,a;GM,x/0. (3.30) 

\ \ 2{ax — pY X J J 

The corresponding solution to Eq. (1.4) in terms of variables S,t,u is given by 



(5,t) = ciS^-^ef-^ +'-S,,= r ^'"".^ + -■ (3.31) 



p ' 2{ax — pY X 
where ci, x G M, x 7^ are arbitrary constants. 
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Case i?3. The last subalgebra hs from the optimal system (3.22) is spanned 

by 

hz=< — +yS— + {xS-xpu + yu)—>, a;,yGM. (3.32) 
We take the two invariants z, w 

p y 

as the new independent and dependent variables reduce Eq. (1.4) to the ODE 

(r — y)zwz — rupw (3.33) 
a^a^ {z^Wzz + 2(1 — Kp)ziUz — npil — Kp)w){zwz + (1 — np)w)'^ ^ 
2 {z'^Wzz + (a + 2(1 — np)zwz + (1 — Kp){a — Kp)wY 

Eq. (3.33) possesses a solution of the type 

w{z) = ciz", ci G M, (3.34) 
where g is a real root of the fifth order algebraic equation 

((r - y)q - rnp){q{q - 1) + (a + 2(1 - Kp)q + (1 - np){a — np))"^ 

2 2 

+^{<l{q - 1) + 2(1 - Kp)q - Kp{l - Kp)){q + (1 - Kp)f = 0. (3.35) 
Respectively the solution to Eq. (1.4) takes in this case the form 

W(Z} = Cl 

gq+i-np^-yqt ci , y , At G M, y 0. (3.36) 

p 

4 Conclusions 

In the previous sections we studied the Sircar-Papanicolaou model (1.3) in the 
case of the nonlinear demand function U{z) = Pz, /3 > 0. The model which 
includes the linear demand function was studied in [14] with numerical methods. 
In this paper we use the methods of Lie group analysis which gives us a general 
point of view on the structure of this equation. Wc found the symmetry algebra 
admitted by the nonlinear PDE (1.4) for r 7^ and by Eq. (3.1) in the case 
r = 0. We present in the both cases the optimal systems of subalgebras. Using 
the optimal systems of subalgebras wc provide the complete set of non-equivalent 
reductions. In most cases we solve the ODEs or present particular solutions 
to them and respectively to Eqs. (1.4) and (3.1). The explicit and parametric 
solutions can be used as benchmarks for numerical methods. 
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